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AIR FORCES AND MOMENTS ON TRIANGULAR AND RELATED WINGS WITH SUBSONIC 
LEADING EDGES OSCILLATING IN SUPERSONIC POTENTIAL FLOW 1 


By Charles E. Watiixs find Julian H. Berman 


SUMMARY 

This analysis treats the airforces and moments in supersonic 
potential j flow on oscillating triangular icings and a series of 
xwepthack and arrow wings with subsonic leading edges and 
supersonic trailing edges. For the icings undergoing sinusoidal 
torsional oscillations simultaneously with vertical translations, 
the linearized velocity potential is derived in the form of a power 
series in terms of a frequency parameter. This method can be 
useful for treatment of similar problems for other plan forms 
and for wings undergoing other sinusoidal motions. For tri- 
angular wings, as many terms of such a series expansion as 
may be desired can be determined; however, the terms after the 
first few become very cumbersome. 

Closed expressions that include the reduced frequency to the 
fifth power, an order which is sufficient for a large class of 
practical applications, are given for the velocity potential and 
for the components of chordwise section force and moment 
coefficients. 

These wings are found to exhibit the possibility of undamped 
torsional oscillations for certain ranges of Mach number and 
locations of the axis oj rotation. The ranges of these parameters 
are delineated for triangular wings. 

INTRODUCTION 

This report is concerned vrith the derivation of expressions 
for the velocity potential and associated forces and moments 
for oscillating triangular wings in supersonic flow. The 
boundary-value problem for the linearized velocity potential 
for an apex-forward triangular wing oscillating in a super- 
sonic main stream may be classified, according to reference 
1. as “purely supersonic” if the leading edges of the triangle 
are outside the Mach cone emanating from the apex of the 
triangle or “mixed supersonic” if the leading edges are inside 
this Mach cone. 

In the purely supersonic case the principle of independence 
holds; that is, the flow on the upper surface of the wing is 
independent of the flow on the lower surface and vice versa. 
Garrick and Rubinow (reference 1) have shown that the 
boundary -value problem for the velocity potential in the 
purely supersonic case can be satisfied by simple distributions 
of sources with local strength proportional to the local pre- 
scribed normal velocity of the wing. Nelson (reference 2) 
has treated the oscillating triangular wing for tliis case to the 
third power of the frequency. 


In the mixed supersonic case the principle of independence 
does not hold. Boundary-value problems for lift-producing 
wrings in this case can be satisfied by distributions of doublets; 
the relation between doublet strength and normal velocity 
of the wing is, however, in general, not simple. The deter- 
mination of this relation requires the solution of an integral 
equation that employs the potential of a time-dependent unit 
doublet as kernel and limits of integration that depend on 
Mach number and wing plan form. 

For treatment of problems that involve boundary condi- 
tions that are independent of time — such as constant angle 
of attack, constant rate of pitching, and so forth — the doublet 
potential, employed as kernel of the integral equation, is 
considered independent of time and in these cases the integral 
equations for triangular wings can be solved by a straight- 
forward process. 

For treatment of problems of oscillating wings, however, 
it is necessary^ to employ, as the kernel of the integral equa- 
tions, a doublet potential that varies harmonically with time 
and in this case the solution of the integral equation, gener- 
ally, becomes very cumbersome. If the doublet potential or 
kernel is expanded in terms of the frequency of oscillation, 
however, use can be made of knowledge of solutions of inte- 
gral equations for problems that are independent of time to 
obtain an expanded form of solution for a wring undergoing 
harmonic oscillations. Such a procedure was demonstrated 
in treatments of rectangular wrings in references 3 and 4. 

The purpose of the present report is to make use of the 
expanded form of the velocity potential to obtain the forces 
and moments, based on the first few terms of this potential, 
for a rigid triangular wring performing vertical and pitching 
sinusoidal oscillations in mixed supersonic flow'. Although 
as many terms of the expanded potential as may be desired 
can be obtained after the first few terms, the process becomes 
very cumbersome. The flow normal to the leading edge is 
subsonic but the flow normal to the trailing edge is considered 
to be supersonic. This latter consideration implies that the 
potential derived for triangular plan forms may be used to 
calculate the aerodynamic forces and moments for other plan 
forms that may be formed with the triangular wring by cut- 
ting the trailing edges so that they lie ahead of the Mach 
cones emanating from their foremost points. 

Other approaches to the solution of the problem of oscillat- 
ing triangular wings have been given by Robinson (reference 
5), Haskind and Falkovich (reference 6), and by Stewartson 


1 Supersedes .VAC A TV 2457, “Air Forces and Moments on Triangular and Related Wings With Subsonic Leading Edges Oscfflattng In Supersonic Potential Flow" by Charles E. 
Watkins, 1SUL 
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(reference 7). In both references 5 and 6 formal solutions 
to the problem were obtained in terms of special systems of 
curvilinear coordinates. Robinson’s solution was given in 
terms of a double summation of trilinear combinations of 
Bessel functions of the first kind with Lam£ functions of the 
first and second kinds. Similarly, the solution of Haskind 
and Falkovicli was given in terms of summations of Bessel 
functions of the first kind combined with elliptic integrals 
of the first and second kinds. In both references 5 and 6 the 
potentials were not reduced to useful forms for calculating 
forces and moments. 

In reference 7 Stewartson makes an interesting though 
specialized use of the Laplace transformation to develop a 
method whereby terms of the velocity potential for triangular 
plan forms, expanded as herein, can be obtained. StewartsoD 
gives formulas that, except for errors presumably' in printing, 
can be used to develop the potential to the second power of 
the frequency, but he omits many details in his derivation. 



SYMBOLS 


disturbance-velocity potential 

*,y,z 

rectangular coordinates attached to wing 
moving in negative z-direction 

Z>v 

rectangular coordinates used to repre- 
sent space location of doublets in 
zy-plane 

z m 

function defining mean ordinates of any 
chordwise section of wing such as 
y~ 2 /j as shown in figure 1 

vi(x,yut) 

vertical velocity at surface of whig along 
chordwise section at y=yi 

Xo 

abscissa of axis of rotation of whig as 
shown in figure 1 

t 

time 

h 

vertical displacement of axis of rotation 

ho 

amplitude of vertical displacement of 
axis of rotation, positive downward 

a 

angle of attack 

do 

amplitude of angular displacement 
about axis of rotation, positive leading 
edge up 

fi,a 

time derivatives of h and a, respectively 

V 

velocity of main stream 

c 

velocity of sound 

M 

free-stream Mach number ( Vic) 

CO 

_ M 2 co 

frequency of oscillations 

w - yp 


k 

reduced frequency (Sm/T 7 ) 

t 

6— tan e 

half apex angle 

®i»m 

represents functions of <3, x, and M 

D n ,F n 

functions used to denote doublet distri- 


bution functions 

Wl 

constants associated with D n depending 
on 0C 


<rj 

2b 


8 

a 


q,N,6,T 

p 

Ap 

P 


M a 




m a 


Po — &C 
F',E' 


constants depending on / 3C 
constants depending on j SC and M 
root chord of wing 
semispan of wing 

slope of ray passing through vertex of 
wing 

dummy variables 
density 

local pressure difference 
section force (total force at any spanwise 
station) 

components of section force coefficients 
section moment (total moment about 
z=z 0 at any spanwise station) 
components of section moment coeffi- 
cients 

total component of damping-moment 
coefficient 

complete elliptic integrals of the first 
and second kinds, respectively, with 
moduli Vl — Po 2 


ANALYSIS 

BOUNDARY-VALUE PROBLEM FOR THE VELOCITY POTENTIAL 

Referred to a rectangular coordinate system moving 
forward at a uniform supersonic speed in the negative 
z-direction (see fig. 1), the differential equation for the 
propagation of small disturbances that must be satisfied by T 
the velocity potential is 



The main governing boundary condition to be satisfied by 
the velocity potential is that the flow be tangent to the 
surface of the wing or 



=w(z, y, t)= 


rr . dZn 

y dx ^ dt 


( 2 ) 


where Z m is the vertical displacement of any point of the wing. 
For the particular case of a wing independently performing 
small sinusoidal pitching oscillations of amplitude cto about 
some spanwise axis z 0 and small sinusoidal vertical trans- 
lations of amplitude ho, the quantity Z n in equation (2) is 

Z n =c i “‘[a t (x — x 0 ) + ho] = a(x — Xa) + h (3) 

(See fig. 1 (b) for sketch showing instantaneous displacement 
of section y=yi.) For convenience, the frequency of oscil- 
lation of both pitching and translatory motion is denoted by 
w. Considering these motions to occur at separate frequen- 
cies would add no difficulties to the derivation. 

Substituting the expression for Z m (equation (3)) into 
equation (2) gives 

w(x,y,t) = Fa+ 6t(x— xo) + K (4) 

Equation (4) implies that the velocity potential may be 
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expressed as the sum of separate effects due to position and 
motion of the wing associated with individual terms of this 
equation, namely 

<£=<£«+<£* 4- <f>* (5) 

DERIVATION OF <f> 

In order to obtain the analytical expression for the 
potential <j>, it is necessary to derive only one of the subsidiary 
potentials appearing in equation (5), say <f> a . The other 
subsidiary potentials 4>* and <&, can then be obtained from 
the derived expression for <f> a by simple comparison. 




(at Plan form u If -plane). 

(b) Section t~yi (ri-plane). 


Flam I.— Sketch iBastntflng coordinate system and the two degrees of freedom a and h. 
272483—54 50 


In order to satisfy the boundary-value problem, for <f>*, a 
convenient procedure is to start with the expanded form of 
the potential of a uniform distribution of doublets. Then, 
for a given power of the frequency of oscillation this poten- 
tial, as will be shown in the following analysis, can be 
modified so that, when integration is made over the appropri- 
ate region, the results satisfy the differential equation (1) to 
the given power of the frequency and satisfy the condition 
of tangential flow exactly. The type of doublet required 
is that with its axis normal to the plane of the wing. The 
potential of such a doublet may be obtained from the poten- 
tial of a source, located in the plane of the wing, by partial 
differentiation with respect to the direction normal to the 
plane of the wing. Similarly, the potential of a distribution 
of the required type of doublets can be obtained from a dis- 
tribution of sources located in the plane of the wing. 

The potential at (x,y,z) due to sources located at points 
(£,j?, 0) in region r (illustrated in fig. 2) of the xy-plane which 
satisfies the differential equation (1) may be written as 



rj8 s FjS 1 


and 

f?= V0r 
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Expanding the integrand of equation (6) into a power series in u, collecting terms with respect to £, and differentiating the 
resulting integral with respect to z gives the expanded form of the potential of a uniform distribution of doublets, namely 


+ • • • +flowS* S,_3 + • ■ ■ ■ • • d'Oinf? Js ^ 3 + • • • ...+ 

r(«"i^+a , t2-ff+ • • • +anmfl 2 “' 3 + • • • 


where 


and 




— \ 2 78-2 


(-!)»-» / 5 V 


a »i 


(7) 


( 8 ) 


(9) 


““ (2m— 2)! 

For convenience in the succeeding discussion and analysis equation (7) may be written hi the following form: 

*-tssWJ ' j? « *] H0) 

An interesting and significant property of equation (10) is that the coefficient of each power of £ satisfies the differential 

/iaiV 

equation (1) and has the form of a source potential with strength proportional to { j^J e~ f3( . This property may bo 
shown by writing the coefficient of £" as follows: 

f 

B L 


“ ( — 
1+ X TXT, 

777 — 2 (2m 


=(!)■<-“{ 


cog 


(M 


( 11 ) 


A more general solution to equation (1) may thus be obtained by introducing properly chosen weight or distribution func- 
tions (denoted by £/?„(£,?;)) into the coefficients of £" in equation (10). Let this solution be denoted by fa; then it 
can be written as 

TT iH " + (12 > 

Examination of equation (12) shows that, at the surface 2=0, the potential fa is determined by the first integral expression 
but that both integral expressions may give rise to normal velocity. In succeeding steps in this analysis it is shown that the 
distribution functions !?,,(£, ij) in equation (12) can be determined so that the first integral expression taken alono will exactly 

satisfj r the boundary condition of tangential flow for <j> a ; that is (j^~) ~^ a - Also, any additional normal velocity that 

arises from the second integral expression can be canceled, to the required order, by consideration of additional doublet solu- 
tions to equation (1). The problem of satisfying the boundary- value problem for the velocity potential <f> a may thus be re- 
duced to that of determining the appropriate distribution functions and additional solutions to equation (1). 

In order to show that the first integral expression in equation (12) can be made to satisfy the boundary condition for 
the coefficient Vaa H i appearing in tins equation is first considered. If the analytical expression for the coefficient a„i 
(equation (8)) is multiplied by x n and summed with respect to n, the result is identically Va. Tliis result may be shown 
as follows: 

(tSz)" 


Va 22 a:"g»i — Va X 


n— 0 


n-0 


n\ 


e -*«*= Vae mx e _<s * = Va 


(13) 


It is significant that this identity holds if only terms in w to any given power are considered. For example, retaining 
only terms including « to the fifth power gives 


Fa 


|^1 — iwx- 


24 




' } m) +x ( is + sjs *- 


24 


M-l 


12 


! )+ 


'(■ 


6 6 


12 


M S-fXm)]-*- 


(14) 
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Next consider the normal velocity at 2=0 associated with 4>i, namely 

T “2 i? g <*'+g (fer=9 © ! "’’“"XJ d.(«,,)£-h— - it *,] (15) 

Examination of this equation and equation (13) reveals that the first integral expression on the right side of equation (15) 
yields 1 a exactly, provided that the distribution functions are determined so that the following integral equation is satisfied: 




ax* 


(16) 


The kernel of this integral equation has the form of a steady-state doublet potential. The problem of determining the 
distribution functions for this case is therefore analogous to determining distribution functions for certain steady-state problems. 
The distribution functions for steady-state problems, at least for those involving conical flow, can be determined by a straight- 
forward process. The main details of this process are given in the appendix, where a method of solving equation (16) for a 
triangular wring is derived and the distribution functions required to derive the velocity potential for this wing to the fifth 
power of oi are given. It is to be noted that the method derived for solving equation (16) for a triangular wing may be 
generalized to apply to various plan forms and to problems of satisfying the boundary conditions for various velocity 
distributions. 

From this point on, the analysis is restricted to the derivation of terms of the expanded potential involving 3 to the fifth 
power. The method of deriving these first few terms is quite general and can be used to obtain as many additional terms of 
tile expanded potential as may be desired. As previously pointed out, however, terms of the potential after the first few 
become very unwieldy. 

If the appropriate distribution functions are known for terms involving 3 to the fifth power, equation (15) may be writ- 
ten as follows: 

w,= Va+Wt (17) 

where 


Va v d 1 3 s 
».■;= hm 




, 01- X 

1 — twr 


Uo I £ 3 '' 


12AP~ 12Ap)^^~12lP d! = dv ( 18 > 


(-* 
V 12 i 


2 1 6 


is the additional vertical velocity arising from the second integral expression in equation (12) involving 3 to the fifth 
power. In order to maintain the boundary condition for <£„, this additional velocity w t must be canceled. As pre- 
viously pointed out this canceling, to a required order, can be achieved by considering other doublet solutions to equa- 
tion (1). For this particular case consider relations & and & similar to (equation (12)) having the following forms: 

’‘-tctSW/ 5,1 *' *>+ SfeSj©’""»-'XX <*»] 

* ,= 75® £ KfX P -«-’> TT d ’+£ feSj) d,] 

The vertical velocity distributions and w t , to the fifth power of 3, arising from these expressions are 
X ^ Jp (1 ds dr, 

I? (^)*XJ[ (1 dt d n 


(19a) 

(19b) 


(20a) 

(20b) 


In [these equations the distribution functions D 0 , D u and so forth can be determined, as discussed subsequently, by the 
method given in the appendix so that Wi+u?* is identically equal in value but opposite in s ign to u’ 2 . When these functions 
are determined, the boundary condition for <j> a is satisfied by Wi+Wi-f w t =Ya, which implies that the potential <j> m to the 
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fifth power of 3 is given by the sum . 

$«=4>i+<£i+<£i (21) 

to this power of 3. 

Expressions for U a (n= 0, 1, 2, 3) and D u (n=Q, 1) are given together with expressions for other distribution 
functions D n (n=0, 1, 2, 3, 4, 5) in the appendix. In regard to the determination of D n and Z5„, it may appear necessary, in 
order to formulate integral equations for these functions, to perform the generally unwieldy integrations of the type 


Urn ~ JJ D K 0-, v )r+W-»dZ dr, 


( 22 ) 


appearing in equation (15). In general, however, the information necessary for the determination of the functions l) H and D H 
can be obtained, as is done in the derivation of the functions D n in the appendix, by examining the values of these integrals 

and their derivatives with regard to the parameter 8=~ at some particular value of 6. 

Returning to equation (21) and introducing into this equation the expression for fa (equation (12)) and the expressions 
for fa and fa (equations (19)), each to the fifth power of 3, gives for the potential <f> a to this power of 3 the following result : 


♦•-^liS'^XXWO 

m[( 


-lux- 


3 s z* . iu t x i , u'x* 


24 120 


■ 3 2 

iu 3 x 

3 4 z* 

i&j'ar’N 

2A P 

2 AP 

4 A f*" 1 

12A/7 


iu-\- aj*a: - 


24 


-)a- 


/ iu 3 u*x iu 1 \ „ iu 1 pj - ] r/to* iu 3 x u*x* , iu 3 x 3 \ 1 

\2il/ 2 "' 2il/* 4A/V" i '24Af 4 " J ^ 2 L\2 2 


st^Ki-**- 


+ 

T~)r~{ 2W ~ 2Ap) B ] + 2M* ^ [( tw + a ' x ‘ ~ l ~T~) 5~2^ s 
2 W ^\_(~2 2 ~ ) 2 M 2 ^ 5 [( V ) 4 W [ (1 ~ tw:c) J + 4M‘^‘ d * dv 

Since it may be shown, as in reference 8, for example, that 


(23) 


equation (23) reduces, at z— 0, to 

<£„= Va£^l — iux 


6 


T“yl?XX^ itrf, - v “" D - ( ’'' , ' > 


) #,Xz,y)+z^iw+w a z — - 


24 

;z» z* 


120 
AMS')** 


24 


12 


+ 


4M sr 12M J 

iu°x\z ? 


)#o (x,y)+; 


i( u * 

iu*x 

\ 2 

" 2 

' 3 1 

iu*x 

2 M* 

2 A/ 4 




U'X IU°X , \ ?=r , 


X 


\4 M 


4 M K 4 M 4 


\2M i ~2M i 4Af*, 
:)^,2 /)+§pA(*,y)] 


0 D, ( 3 ‘ ,V) (4 M* 4 Af 1 ) Dl 12 A. p D &> y )+ 


(24) 


which, after the expressions for the distribution functions given in the appendix are substituted and the terms are 
regrouped, may be written in the following simple form: 


fa = Va^C i x t —y s [A 0 — iux v 1 —u s (<r i x i -\-<r t p i y t )-\-iu 3 (<r B x 3 -\-(r ll fl t y I x)-{-u\<r t x t + tr^^x* + om/Sy) — iu s (<r uX*+ 
<ruP*y*x*+<TtiP*y*x)] 


( 25 ) 
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The quantities A 0 , <r u and so forth in this equation are part 

0 R 

of a group of quantities A t and which are 

defined in the appendix; these quantities are suitable for 
writing the potential and expressions subsequently derived 
for forces and moments in simple form. The quantities A t . 
P } , Qj, and Rj are functions only of the product /S C (ratio of 
tangent of the half apex angle of the triangle to the Mach 
angle) and are tabulated for some values of jS C in table I 
and are shown plotted in figures 3 and 4, respectively. The 
quantities <rj, which are functions of Mach number Af and 
the product ( 3 ( 7 , may be evaluated for particular values of 
M from the values of Pj, Qj, and Rj in table I or from the 
plots in figure 4. 

The quantities A 0 and Ai are the same, as should be 
expected, as the parameters associated, respectively, with 
constant angle of attack and constant pitching of triangular 
wings of references 9 and 10. 

Expressions for the potentials fa and 4>k can be obtained 
by the method discussed for obtaining 4> a , or they can be 
obtained to the fifth order of u by comparison and synthesis 
from equation (23). After simplification these expressions 
are 

<£*= d jv'CV — '/ — i ufo-iX 2 — ri-aiST/ 2 )— mVt**— + 

cru0 y* 2 + ott^V) +“' Vis * 1 + viofS : , y l r‘+ 

Vm/ 3 4 »/ 4 t )] | ( 2 6 ) 

<l>i=K\C*x t —y t [ri, 0 — fwro-i— u 2 (o'*^ I + v 4 ^ 1 y*)+ 

iu 3 (tr s x s + tr,fi 1 y 1 x) + w\a- i x i + o-iofSV * J + onjSV)] ( 2 ~) 



At BC= 1 or C=\i which is the condition at which the Mach 
P 

lines from the apex of the triangle coincide with the leading 
edges of the triangle, equations (23), (26), and (27) reduce, 
respectively, to 

K5+7AW- 

(5-2AP)M+2^3? K7+3AP)x*- 

(7— 2lfW*l + 3 = 3 ^* [(83lP+414AP+63)r 4 + 

42(2A/ 4 -9AP -3)j9V x 2 +(8lf 4 -36 AP+ 63)/3V] - 

349 II 3 J 4 [<73AP+638AP+297)r s +2(46AP- 

253AP— 297)j3yr*+3(8AP— 44AP+99)/3Vr]| 

(28) 




iK3 


■■■■■ 




l«IB»’l 


HHiiiiiiiiiiii 
■■■■■■■■■■■■■■■■ 


limn 




I.IBBBBBBBBBI 


■■■■■■■■■■■■■■■■■■■■I 


9IIIIUHI1I 


IBBBBBflBBBBBBBBBBBBBBBl 


^MSSSnSaSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS^H 


IlKSSl !■■■■■»■■■■■■■■■■■■■■■■«■■■■■■■■ 
BBKl iBBBHBflBBBBBaBl iBBaaaBBBaaBBB 

IDIssiisiiiiiiiiisiiiiiiiiiiiiiiiiiiiiiil 

BVBflBMBBBBBBaBBBBBBBBBBaaaBBaaBaBBBBB 
BaBBaaaaaaaaiBBaaBBaaaaaBaHBaiMBaaa 
: aBMBBaaBaaaBBBBBBBBaaBBaaBaaaaBaaaaaa 

:' : -i <' ibbkbbblibbbbbbbbbbbbbbbbbbbbbbbbbbbbbb 
aBa.iBBaaaaaBBaaaaaBBBaBBBaiaaaaaaaBBB 
^■aaavBBBLiBaBBBaBBaaaBflaaaBBBaaaaBBaaaB 
laaaaaaBBaaa'aBaaaaaaaaaaaaaBBBBaaaaaaaaaa 
jBBBBBBBVBaaBaaBBBBSBaaaaaBBBBBaaBBaaBBBa 
BaaaaBaBBt'aBaaaBaaaaaBBBBaaaaaaBBBBBaB 

BBBBBaaaaSvaB BBSBBBBaSaaaBaBaaBBBaa 
■angBBaBUBaBiBBaaaaaBBBBaiaBaaaBBBBBBaa 

■aBCBBBaaa^BBABflaaaaBaBBBaBBBBBBBBBBBBB 

I IHBflBk?aaaBaBXBai.iBBBBaBBaBBBBBBaBaBBBaBaa 
KiBaaBgBaBaaaaBBMaaBBaaaBBBBBBBBBBBBBBB 
iBK&a&aBfiaaaaBaBaKBBaaaBaaaaBBBaaBaBBaBaB 
liigaaBBBfiaBBaBKaa^aaaaBaaaiaBBaaaaBBaaBB 
IfiH^aBBaaKaBaBBKBatiiaBaaBaBBBaaBBaBBaaaB 
^■aaaaBBBtiaBBaaKaBSaBaaaaBBBBBBaaBBBBaa 

■ "~BBB?aBB8aaB£SiaaaBSBaBaBBBBBBBBBflBBBfl 
^^■BaaacaaBBaaafliisiaBaasasaaaBBaBBflBflaBa 
^■iaSSBBKi!KBaBBBBBBkS!IBaB£!<SBBBBBBBBBBflg 
MaaaaaBBgaaBBBBBaaaaBgaBBBggjiBaaaBB— I 


=sp 

Ibbbb 


mi 


: ?BBB BBBBBBBBBg»«BBggBBBBBBBBB»B5 

— MBggPBBBBBBBBaBC BBBBaggBBgga 

fiBBBBBBBkBbi — MB— Bl 


Z 3 


.8 3 





870 


REPORT 1009- NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


(<#>a)3C-1 


, 4«Va: a — p*y 
3|3x 




-UV)- 



icox w s 
~3 


[(5+7il/ J ^ 2 - 


2^ [(llM s +7)^-(2M , +7)/3V^]+ 

l2(13Zf 1 +27)x*+3(M*-21)/3yz t - 

(4tf-9)^y < ]+ 4 158 ^ 7< [(457M 4 +2046il/ ! + 
297)2 5 +2(118il/ 4 — 1221i/ 2 — 207)iS s 2/*3 8 — 
(168A/ 4 — 396M 2 — 297)/3 4 y 4 i] J (29) 



(5-2M 2 )«+ 21 g 7 , [(7+3A/ 1 )®*- 

(7-2il/ i )/3y 2 ]+ 37 g|^J(83M 4 +414il/ 2 +63)x 4 + 

42(2A/ 4 — 9M 2 — 3)f} , y , x i -f- (82/ 4 — 362/ J + 63)0 4 y 4 ] J 

(30) 


Equations (28) , (29) , and (30) can also be obtained by going 
to tbe limit /3C=1 in the expanded potential for triangular 
wings with supersonic leading edges (see, for example, refer- 
ence 2 for results to the third power of «) or by integrating 



(0) At. 

Finos* 3.— Continued. 


Cd) At. 

Fioor* 3.— Continued. 
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the expanded potential of a unit source over the region oc- 
cupied by the wing with 0(7=1. Equations (28), (29), and 
(30) therefore serve as a check on the results given in equa- 
tions (25), (26), and (27), respectively. 

For values of the product 0(7 such that 0 1 (7 i -Cl, equa- 
tions (25), (26), and (27) reduce, respectively, to 

4>*=V ayC'x'-y 1 (31) 

<f>i=a^C 2 x 2 —y t (_x — x 0 ) (32) 

(33) 

These expressions are the counterpart of the potential associ- 
ated with triangular wings of vanishingly small aspect ratio 
in steady flow (reference 11) and may thus be regarded as 
the potential associated with a small-aspect-ratio triangle 
oscillating in pitch and vertical translation in either sub- 
sonic or supersonic flow, provided that 3 remains finite. The 
condition 0(7*3; 1 is obviously satisfied for all values of C 
when M=l. However, as A/-*l, the value of 3 becomes 
infinite and the expanded potential, as treated herein, be- 
comes meaningless. 

FORCES AND MOMENTS 

As pointed out in the introduction, the velocity potential 
for the triangular wing can be used to calculate the aero- 
dynamic forces and moments for other plan forms that can 
be formed from the triangular wing by cutting the trailing 



le) Ai. 

Fir, i; as 3— Conclude*!. 


edges so that they lie ahead of the Mach cones emanating 
from their foremost points. Sketches of different plan forms 
thus obtained are shown in figure 5. 

The force and moment coefficients desirable for most flutter 
calculations are those that yield the spanwise variation in 
these quantities or ehordwise force and moment coefficients. 
These coefficients are obtained by integrating the pressure 
difference along any chord for the forces and the pressure 
difference multiplied by a moment arm for the moments. A 
convenient procedure in deriving these quantities is to intro- 


duce the reduced frequency parameter ^—k and to employ 


the variables x, y, and x 0 in a new sense as nondimensional 
quantities obtained by dividing the old variables by the 
maximum chord 26 of the wing. 

The pressure difference between the upper and lower sur- 
faces of the wing is 




(»> Pi- 

FE'iURE L — Variation of the quantities P Jr and /?* as functions of pt 
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Figure 4— Continued. 
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Figure 4*— Continued. 
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The force, positive downward, at any section of any of the 
plan forms shown in figure 5 may be expressed as 

p - 2} ],'> d '— teb i‘,',c(jc s? + W d * 

= -2 p fUT 1 +2ik[ 11 4>dx\ (35) 

( J Vic Jtlc ) 


where Xi has the following values with respect to the different 
configurations shown in figure 5 : 


In plan form A 
In plan form B 

Xl=l 

(36a) 


/!= 1— — 
m 

(36b) 

In plan form C 

m 

(36c) 

In plan form D 



a*i=l+— 

m 

* 1= 2 fe 

for O^ygm(^-l) | 

,orm GtC 0 <!,S 26 ( 

(36d) 


After the expression for <f>, given by the sum of equations 
(25), (26), and (27), is substituted into equation (35) and 
the integration is performed, the results may be reduced to 
the form 

P=-4p6^V“‘[^(A 1 + ti a ) + a fl (L 1 +ii 4 )] (37) 



Po 

(1) R,. 

Fkiure 4.— Concluded. 


where 


i 1 = VC i *. , -y 1 [|5(2MV 1 - i 3 J rlo)-^^(4ilP ( r 5 -^V ! )-^^(8M ! CV 8 + < r ! -4 J 8*CV 4 )]+ 

7 / cosh- 1 ^ ( ff|+4 p*C*<rS] (38a) 

L i =\C*x^—y i (3iTW ff «+<r 1 )+— ^ ( f - (5 A/V,-^V t )+ 

(15 M 2 f} 2 C*ff n + 2 a j + 5 ^ , C*<r e ) J (38b) 


U-J&xW \jp+J^ (M 2 i8V 1 -3MV,+3M*/JV J - j 8«A 1 )-^j^ (3M^ J < r 4 +il/ 2 <r 1 + 3ilf ! /3 2 C^l,- J 3M 1 )+ 
16i !p-" - " {5 M«4r,-W«r.-SJi , |8 , a 1 ,+p« gT )+ 16 ^^ V (15 **,.+ AP^-oAPpK’^,- 1 5 MW<r le - 
^Vr5^ < t7 2 <r 8 )+^j||^(15M 4 i9 s (7Vi 1 +2MV 5 +5M s i3 ! C ,i t r s '-15M ! /3 4 C , V lT -2i3VT+5/3 4 <7 s <r*)J-2a; 0 L 1 (38c) 
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L^C'xf-y* (p'A a + 2(3'A 1 -2APa{)+ 2 —t X ' Z (4M^ s -AP^ t -4AP^a r +^a 2 )+-j ^ (8Att7V,+ JfV- 

4APf}*C** i +8ZP(i i C* t r i -(} i r i -4fi*C*Aa-± ^^ &M t <r a -M i fa t -6Wp** u +p* 1 j- - ^ X ££ (MMVir a + 
4M-<r i -24AP^C I ^-96AP^C i <T lt -4^ a - lli +24^ i C i l r l ^-^^^(Z2M^ t C i ffli +2AP^+4AP^C i <r, 0 - 
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0 1 <r 1 +4^C 2 A t )+=~ i f (23Po- 9 +42/ i jS i (7 1 o- 1 ij+16 2/ 2 (3 W t i— 2/SVii— 4/S 4 C !i £r t ,— 16/WV ir )l-2jr 0 £i (38d) 

In a similar manner, the moment (positive leading edge up) about the axis x=x Q is 


M a =—4b 2 f 1 (x—Xn)Apdx 
Jrlc 

= -4pFW (J/t+tJjf^+a^i/a+iAf*)] 


Af t =yO V-2/ I [|^(2AP ffl - J 3^ a )+^(AP ( r 1 +/3Uo)+^^ 0^,-4 tf V,_ 
5 ^ t ^+APa i + 10 APp i C i <r t )-^^^( 2 p i cr i + 5 ^C 1 <rt + 2 .\Pa s + 5 AP^C a aSy 2 x o L l 


M^CW^y* [^p+ ^f - Jl> 03V,-3 APaJ--^£g£ (fr l +APr s +4APf(Pa i )- 1& ^ W (pa 5 -5AP* t )+ 

(0V.-6/* <CV.+3fV # + 18W^,o)+ 2 ^ lif< (/JV.+JWr.+AfV^WCW 

83/WV„)]+y 2 C osh-^[^-^ (iV/ J vj+4il/ J ^ 2 C 2 ffi+/5Vi)+ 

(l/=v s +2.I/ 2 |3 J C I < r l o+8^/S‘CV I1 +^V j +2i3 < C s v 5 )]-2r Q i 1 (40b) 


M i =fC'x l '-y i OTV, - 33fV,+ 3 M 2 £V S - 0^ t )- ^ (W<r 1 +AfV a +43/^*(7 , <r*- J 9U 1 - 

1 ft \/4Z.*- r fi J. 3 -ij* 

3f 2 dV 2 +43/ 2 g 4 6’M 3 ) — 3 i3 t - (3/^V 5 -53/V s -^V 7 +S3/^V II )+^|^^(xIf^V i -63/ J /3 4 C'V s +3/V s + 

l 8 M t prPa ia -p i vr-$P*C 1 v i -APf} 1 a ls -l 8 APl 3 V 1 r l J+ 2 — ^ lV * (AP(} i <T s + 2 APp i C t <T t +AP<r t + 2 i\Pp t O‘<T l < > + 

S3/* ( 3^V„-^V T +2^ 8 C'V 8 -3^V lJ -2lf 2 ^^v I! -83^ 5 CV I7 )]+y 1 cosh- 1 ^ (3^V x +3/V*+ 

4A/*/3*C ri <r«~/3U 1 -34 1 j3 J <r 2 +43i 2 i3 4 (7 1 ^l 1 )+^S^(3f , ^ 2 v s +23i , j5 4 CV s +3/Vo4-23/^ 1 C , <r l0 +83/ 4 ^v, l - 

pc 

d V r + 2^ <5 <7 1 fl-s— 3/*/3*<ri S — 2 APfMPcr lt — 8 3f s jS*C Vu) J -2x 0 (3f I +L a +2r a L l ) 


(40c) 
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[||£ (^- 221 ^+ 2 ^ 0 - 

3 ^ 16 ^ a?1< OTV,- 

OTV,~ 

Sil^^+MVs+lOil/^^-jSVj-S^C 0 ^- 

WV+ lOil/^CTVsH (2 W<r,+ 

5lf 2 (3 4 C i < T4+ 2M*<r s + 5.M 4 j3 2 £7V e - 2£ 4 <r 2 + B^CPA,- 

2M 2 !8V 7 +5itP)3^ gg )+ - ^fo 3 - 8 W»r 

7lf 2 i3 4 (? ffw +2l/V 12 +28M 4 /3 2 C i o- I8 - i 8 4 (ri 5 + 

7 pCPcu- M z pa ls - 2 S^iS^cr,,)- 

(4M’£ V 9 + 7 AT/3 4 (7Vo-- 3 5 M 2 iS°g 4 g- 11 + 

4MV ls +7A/ 4 /5*C 2 ( r 1 3+70M 4 )3 4 C' 4 ffu -4j3V I6 - 

7/3«C , *<r 14 +35(S 8 C 4 <r I7 -4A/ 2 )3 s (ri s-72li 2 /3 4 C 2 <r lt - 

70AW?V M )-^|^ (8 W<r,+ 14M 2 /9 4 (? 1 <r 1 o+ 

353f 2 /3 8 C ,4 <rn+8M 4 (r 1 s+14M 4 /3 J C 2 <r 1 j+35M 4 /3 4 C 4 (r 1 4— 

8j3V ls — 1 4£ 8 (7V M — 3 5^ 8 (7 4 ff i7 — 8 AP/SV^ — 

1 4 A/ 5 £ 4 C7V la - 3 5 jtfVCVso)]- 2 z 0 ( Af j +Z 4 + 2 z 0 Z 2 ) 

(40d) 

In equation (37), for example, the quantity (Z, +iZj) is the 
lift-force coefficient associated with vertical motion of the 
wing. The real part Li is in phase with the vertical position 
of the wing and the complex part Z, is 90° out of phase with 
this position. Similar definitions apply to the lift coefficient 
(U+iLt) associated with pitching motion and to the moment 
coefficients (M\ and (Mj+tM 4 ). The imaginary or 

out-of-phase terms determine the aerodynamic damping 
associated with different wing motions. 

Although the expressions for the components of lift-force 
and moment coefficients in equations (38) and (40), respec- 
tively, are lengthy, they may be quite easilj' evaluated with 
the aid of table I or the graphs in figures 3 and 4. 




(a) Plan form A. (o) Plan form C. 

(b) Plan form B. (d) Plan form D. 

Fiodse S . — Sketches Illustrating different plan forms for which the forco equations (38) and 
moment equations ( 10 ) apply. 

DISCUSSION 
SAMPLE CALCULATIONS 

In order to give some indication as to the general nature of 
the spanwise distribution of the different components of lift 
and moment coefficients, equations (38) and (40) have been 
evaluated at different spanwise positions y for plan form A 
(fig. 5(a)) for the following set of conditions: £(7=0.5, 
(7=1.0, z 0 =0.6, M=V5/4, and £=0.1. These sample re- 
sults are plotted as functions of spanwise position in figure G. 
The spanwise variations of the different components of lift 
force are shown in figure 6(a) and the corresponding varia- 
tions of moment coefficients, in figure 6(b). 

In figure 6 note that, for the particular set of conditions 
for this example, the maximum numerical values of the com- 
ponents of moment coefficients M lt Aft, and M 3 are positive 
and are near the tips of the wing, whereas the maximum 
numerical value of the component Af* is negative but is also 
near the tips. It may also be noted that the integrated 
(in spanwise direction) values of the components of moment 
coefficient, or components of total moment coefficient, w’ould 
in each case have the same sign as the maximum value of the 
corresponding component of section moment coefficient. 
This result is not necessarily true in general, because changing 
some of the parameters involved in the evaluation of the span- 
wise distribution of some components of both force and 
moment coefficients may change the distributions signifi- 
cantly from those shown in figure 6. 

The fact that the total component of moment coefficient 
Mi is negative in the 1 example just discussed shows that, for 
the conditions of the example, this term would not contribute 
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iaj Lift force. 

BVtUK* 6.— Span wise distribution of components of lift force and moment coefficients for 
n-0.8. tfC-0.5;C-U0;t-0J; Af-VS/i 

to the aerodynamic damping but, on the contrary, would act 
as a source of energy for the oscillating system. This cir- 
cumstance is significant since it leads to the possibility of the 
single-degree-of-freedom torsional instability discussed in the 
following section. 

UNDAMPED TOESJONAL OSCULATIONS 

The wing plan forms discussed herein, like two-dimensional 
and rectangular wings, exhibit the possibility of undamped 
torsional oscillations for certain ranges of Mach number M 
and location of axis of rotation x 0 . This fact is borne out, 
as indicated in the preceding paragraph, by considering the 
integrated (spanwise) value of the component of damping 
moment M t associated with pitching or torsional motions. 





The main results of this phenomenon can be obtained by 
considering very slow oscillations so that only terms in 
equation (40d) for 3/* involving the reduced frequency k 
to the order l/k need be retained. In this case, 

v iy *+^-° 1)^0 

( 41 ) 
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For the triangular plan form (fig. 5 (a)) the integrated value 
or total component of moment coefficient is 

{ -(2M ! +l)io+3(2F-l)i 1 + 

4z 0 [Ao-{2AP- l)^L 1 ]+42 0 s (ikf 2 — 1)4 0 } (42) 

In general, the condition of torsional stability or instability 
depends on the sign and magnitude of M t in equation (42). 
Positive values of M t indicate stable conditions and negative 
values indicate the possibility of torsional instability. 
Between the stable and unstable conditions — that is, when 
M« vanishes — a borderline state of unstable equilibrium 
separatingdamped and undamped torsional oscillations exists. 

The ranges of values of Mach number M and location of 
axis of rotation x 0 for which vanishes for some selected 
values of C— tan e are shown plotted in figure 7. The 
regions inside the curve in this figure indicate instability. 
The dashed curve, on which some of the solid curves ter- 
minate, represents the locus of values of M and x 0 for which 


M< vanishes for the whole class of triangular wings with 
supersonic leading edges; that is, for triangular wings where 
pC^ 1 (see fig. 5 of reference 2). It will be noted that M t 

3 

vanishes for values of x 0 ahead of the root ^-chord position. 


It will also be noted that, as the vertex angle e=tan -t C 
decreases to 30°, the range of values of Mach number for 
which M4 vanishes decreases sharply. 

In conclusion, investigation of equation (42) shows that, 
for a given value of the reduced frequency k, Mach number 
M, and location of the axis of rotation x 0 , the magnitude 
of the damping coefficient Mi generally decreases as C de- 
creases and, consequently, torsional instability is less likely 
to occur with slender triangles than with wider triangles. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., June 19 , 1951. 



Fiouri 7.— Curves showing ranges of Mach number M and ails of rotation zt for which the aerodynamlo torsional damping moment vanishes for some seleclod values of angle*. 
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APPENDIX 


METHOD FOR DETERMINING DISTRIBUTION FUNCTIONS 


In tliis appendix a method of obtaining the distribution functions is developed in detail. Expressions for the distribution 
functions required to drive the velocity potential to the fifth power of the frequency of oscillations and a set of functions 
useful in 'writing the expression for the potential in concise form is given. 

As indicated in references 9 and 10, where the triangular wing is treated for constant angle of attack and for constant 
rolling and pitching motions, a convenient form of the velocity potential corresponding to a distribution of vertical velocity 
proportional to r n is 


. _1 ^ f*» i? +t dZ 

r d 2 J-c V(x— 


(Al) 


where c is the slope of a ray passing through the vertex of the wing (ij=<r£), £1 is the least value of £ that causes the 
denominator R in the integrand to vanish, and F.(tr) is the distribution function that is to be de- 

termined so that 




(A2) 


In equation (Al), the integrand is noted to be singular at the limit £=£i. A form of the integrand which avoids this 
difficulty is obtained by making the following change of variables: 




e=l 

X 


(A3) 


r=cosh -1 

where r is the new variable of integration. With these substitutions equation (Al) becomes 

*•-; HI '* (<- N cosh r) ‘ + ' ir 

and tlie corresponding expression for > which is an Integral equation for is 

“S { ' +: 2) J - c ( 1 V A'*(^-co S h r)” +, co a hrJ.+ 

p*z 2 (n + 2) j_ c Ft j n cosh (jv+ n cosl1 r ) cosil T dr \ 


(A4) 


= x * (A5) 

In this equation the value of w A obtained by performing the indicated integrations and then going to the limit 2=0 is the same 
as the value that would be obtained by first going to the limit and then perfo rming the integration but neglecting singularities, 
pointed out subsequently, that arise when the value of c approaches the value of 0. Making use of this fact reduces the 

S79 
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calculations involved and simplifies the integral equation for F n . Thus passing to the limit gives for equation (A5) 




p%n + 2)C c F n (a)da (*•»**& Ar ,/ q , V+ 1 , , 


=x n 


(A6) 


where A r a is the value of N at 2=0. The requirement that the normal velocity be proportional to x n and independent 
of y implies that all derivatives of equation (A6) with respect to y or 8 must vanish or that F„(a) must be so determined 
that the final value of the integral in this equation be independent of 8. The requirement that all derivatives of equation 
(A6) with respect to 9 v anis h leads to other equations for F a and, after n+1 such differentiations, the equations 
acquire forms for which solutions are known. The value of F n can then be determined from these known solutions 

• d K w • 

by evaluating each of the derivatives a * ^3" arbitrary value of 8 in the range — as follows: 

The ifth derivative of w, with respect to 9 (equation (A6)) gives the following integral equation for F n \ 


d K vj n 


/9*(n+2) f ° F n {a)da /3* V 

x J. c (l-m U 


de* 


fa) (.q-No cosh r)»+ lCO ^$(r)dr=0 


(A7) 


where 


#(t)= 


(n+1)! , 

( /5<r+cosh r \ 

* Kin+l)l , 

( |8(r+COsh r \ 

(n+l-Z)! 1 

cosh r—q) 

1 (n+2— Z)! 1 

V-iVo cosh r—q) 


K(K— 1) fa+ 1)1 / ' j3<r+cosh t V 
(n+3 — K)\ \2V 0 cosh r— q) 


, tt-H 


N' 0 ~ 2j T‘ - - +( — 1) K KIN 0 ~ K: 

When K—n-\-l, the expression for <t>(r) may be recognized as being a binomial expression, namely 

_ ( n+l)ltf«rlVo+g)" +I 
_ (-iVo) n+1 (g-iVo cosh r) n+l 

(n+l)lx n+1 

~(-NJ n+ Kq-N 0 cosh r)“ +I 

Thus the integral equation corresponding to the n+l derivative of w n is 


(AS) 


d a+h w n (— l) tt+1 (7i+2) 
dd n+1 * 


£ 


FnW) da 


iVo B+2 (l-^*cr 2 )” + 2' 


psiDhr 1 

,n+l J 0 




A r «Vi-pv ! cosh t dr 


(— l)» +1 (n+2) Ia»Vl— P 2 # 2 re F n (a) 


r c F n \ 

J_c(0— < 


.^71 -f 3 


da -— 0 


(A9) 


Further differentiation of w n leads to other equations involving integrals similar to that in equation (A9) and, as will be seen 
subsequently, is not necessary for the determination of F„. The singularity at a=9 in equation (A9) is a result of going to 
the limit z=0 after equation (A5) and, as previously implied, is to be ignored. 

Consider the following equation similar to equation (A9) : 


X 


c Ma)da 
-c (9 — a) n+3_ 


=0 


(A10) 


It is known by analogy with problems in incompressible flow and may be shown by direct substitution and reduction that 
this equation is satisfied for any value of n (n=0, 1, 2, . . .) by the function 


/ 0 (<r)=VCW 


(All) 


d K w. 


This function also satisfies equation (A9) for any value of n and satisfies equation (A7) for all derivatives beyond 
the nth derivative. A more general form of solution to equation (AlO) may be shown to be 




(A12) 



TRIANGULAR AND RELATED WINGS WITH SUBSONIC LEADING EDGES OSCILLATING IN SUPERSONIC PLOW 881 


where in is an integer and m^n. The validity of this solution follows, since 


and 


77) | f} m ~ T 


(A13) 


a- m \C 2 — a 2 » ml d m ~ T -s'C 1 — a 1 

(0— <r) B+3 “ ^ (m— r)!rl (fl— o-)”- r+3 


(A14) 


Each term of the summation in equation (A14) is found to have the form of the integrand in equation (AlO), from which 
the function (A12) is concluded to satisfy equation (AlO) and equation (A7) for all values of K^>n. 

From the foregoing discussion and consideration the distribution function F n {a) may be uniquely determined in terms 
of expression (A12) as follows: 

Consider the expression ■ 

= . . . +5n/» (A15) 


where the coefficients B Q , B u . . . B n are constants that are to be determined. 

Each term on the right of equation (A15) is noted to satisfy equation (A7) for all derivatives beyond the nth and a 
total of 7i+l parameters is to be determined. If expression (Al 5) is introduced into equations (A6) and (A7) and the indi- 
cated integrations in equation (A6) and in equation (A7) for K= 1, 2, ... n are performed, n+1 linear equations are 
obtained in B a , B u . . . B n from which these constants may be determined. 

The integrations with respect to a in equations (A6) and (A7) are in general difficult and tedious to perform ; 
however, as previously pointed out it is only necessary to perform the integrations for some particular value of 6 in the 
range The integrals have their simplest form when 0=0 and the integrations can be made for this value of 8 

by reductions and use of formulas in reference 12. 

The functions D n and D a may aIso_be determined by the method discussed for determining the functions D n . The 
conditions to be satisfied by D n and D x , however, are determined by the value of w 2 , equation (18), and its derivatives 
with regard to 0 at 0=0. _ 

The distribution functions D n (n= 0, 1, 2, 3, 4, 5), D„ (n— 0, 1, 2, 3), and D K (n= 0, 1) calculated by the foregoing 
method are as follows: 


D 0 =A 0 \C 2 ?-v l 

Di=Ai VC 2 £ 2 — 17 2 

d 2 =(a 2 +a s 

D 3 =(a,+A s 
d<=(a s +a 7 ^+a s 
D s =(a 3 +A 10 £++A n ^0 2 e-v 


D 0 =\ 


e ) 

\ \C 2 ?-y 

2 

DH 

^A 2 -\-A 3 

r ) 

l Jc'e-v 

2 

v 2 =\ 

(a^+a, 

pv 

? 


\c*e 

B.H 


e 


v 0 <*?■ 

^o=| 

^■Ao+-Ai 

@ 2 y 2 

e 





PV 

e 


\c*e- 


where, with the notations 


ttt, (n + 2)/3 i! C+ c a\C 2 -a 2 , ,. £>* fmar<£ / q V +l u j 

Ir: - — <x-w* d " la sfJ, (s;- C08h v coskTir 


02 //T2 -2 Xp Crnsh-t g 

^ lim ^ ^ » (q—N 0 cosh r)^ 1 [2 ff -(n+4)iV 0 cosh r] sinh 2 r dr 

XX J-C 1 — jSV 9-iO ov y J 0 

and 

1 fi ? r%VCWVl-^V d<r lim ^ N 0 \q —N 0 cosh r)» +1 [4g-(»+6)lVo cosh r] sinh 1 r di 

XX J -c g—iO Oa r J 0 


(Note that the above expressions will vanish identically unless p and m are either both odd or both even.) 
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the .4’s are given by the following equations: 


A 

°~Wl 


A x - 


0,0 


1,0 


AiW^ 0 ~hAtW%t= 1 1 

A s Wl o +A t Wl a =0) 

AiW < s 1 o~\~AiWg t2 = 1 ) 

AiWU+AiWl^o) 
AiWio+ArWU+AsWl^l' 
AiWU+AiWU+AtWl^O j 

A e WU+A 7 WU+A s WU=ol 
A,Wl^A w Wl,+A n WU= 1> 
^L 9 W^ o +^ lo TF! lJ +^ 11 Tn4=0 
A i Wt l0 -\-A w Wi lt +A n Wt,i— 0; 

AoWU+^Wi^AM.o) 

AJVU+AWlrAtWlJ) 


A a Wl 0 +A s Wl*= 

AtWi'O+A 


AtWl^AWU) 

A'Wl^AiWlo] 


^*^2,o+-3 s ^. !t +2.T : F2.*=^,Wi. (1 +^ s Fs. s 

A t Wl, a +A,Wl. a +Z'Wl,<=A,Wl 0 +A s Wl a 

A t WU+AtWU+A a WU=A a Wt' a +A a WU } 


wu) 

^7 yy i,o-t-^i* yrs,t-t-si a wi it =ix l w 5.0-r^-sl^lt ) 

^WU+A^WltAA.Wl^A^Wt^A.WU) 


^ 7 Tno+^ 8 Tn»+^^ T n.4=^Tn 0 +^jri. a 
A 1 Wl a TA,WU+A t WU=AjWl 0 +A t Wh 
A 7 Wi, 


AoWU+l.wit+ltWi^^ Wu+IM.t+AtWU 
^oWU+Ijvu+l.wu^^Wto+A^Wu+^WU 
loni 0 +ljn 2 +zjn 4 =-^tf3.o+^Wi 0 +^ 1 H 7 i 




3,* 


I 8 nio+^ 4 Tii 3 +IjTT 1 =-^^?.o+-Wio+A,rri. f 
A s Wl, a +I t Wi f +2 S Wl .= Wu+AtWlo^UWl 
Zwu+i<wu+Xwu= ffj.o+^Fj.o+^irj,, 

where, with pC=p 0 , 

WU=E' 

Po ^+(I-2po 2 )^ 

Hl ’°~ (1-Po 2 ) 

™ (5po s — 3po 4 )i ?, +(2 — 10p 0 2 +Gpo 4 )£’ / 
n ' 3 -°~ 2(1 -poT 

TTTn 2p 0 ^-(po J +Po < )S / 

n ^~ 20 2 (l-p ( > 2 ) 2 

Wi. a =% WU 


WU= 


C 1 

(5p 0 4 -3po 2 )F / +(6-10p 0 i +2p 0 < )E / 


(1-Po 2 ) 2 


„ ro (27po*-31po 4 +12po e )/ , '+(6-55po 2 +65po‘-24po 8 )S' 

" *•“ * 6(i-^ “ T 

tt^o . (9po 4 -po 6 )i ? '-(3po 2 +7po 4 -2po' , )S' 

M_ “ " 'mi-PoV * " ‘ ' " 


Wlo=%Wlt 

T7 7a (2p 0 6 +9po 4 -3p 0 2 )F'+(6-15po 2 +5p () 4 -4po 5 )S' 

‘ 3I_ (1-Po 2 ) 8 

rrm (56p 0 , -92po 4 +72p 0 e -20po 8 )^ , +(8-117po 2 +202p 0 4 -149p 0 s +40po 8 )E' 

l ' Uo “ 8(1 -po 2 ) 4 

ttto (17po 4 -2p a 8 +po 8 )i ? ’ , -(4po 2 +18po 4 -8p 0 8 +2po 8 )g / 

8^*(1 — Po 2 ) 4 

it t o (8po 9 4■8po 8 )F , — (po 4 +14po s ^-po 8 )•^ , — 

80 4 (1 — po 2 ) 4 ~ 

wi l0 =^wu 

T (epo^SfW-SOpoHepo^'-^-Sgpo’+epo^SQpo'+^po 8 )#' 
w *- a= 
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_ (3p a 4 -6p 0 4 +51p 0 8 ).F'-(6p a *-24p 0 4 +54po ! -R2p 0 8 )£;' 

' 2^1-po*) 4 


w\: a =^wu 

wu^wu 


_ (168p 0 8 — 276p 0 5 -|-216po i — 60po*)i ?v +(24p a 8 — 351p 0 6 +606p # 4 — 447po*+120).E' f 

(1-Po 2 ) 4 

Tr . 0 _t400po 2 -833pa 4 +994p 0 a -553po 8 +120po 1< ’)i ?T '+(40-859po 8 +1910pa 4 -2115pa s -)-1136po 8 -240po to )j£' 

40(1 — po*) 4 

_.. C115po 4 -2p 0 8 +19pa 8 -4p a 1 °)F'-(20p 0 8 +151p a 4 -74p a 8 +39p i> 8 -8p a 1,> )g / 

*'* 400X1- Po 2 ) 4 

ir „ „_ (55po l> +74p o 8 -p C) ta )F f -(5p o 4 +l08po 8 +17p 0 8 -2p 0 ll) )g / 

40/3X1 -po 8 ) 5 

M «.0 fl% ** 


(8 Po 10 — 35po*+ 112po 8 +49 Po 4 — 6po 2 )i r, +(12 — 47 p 0 s — 22 p 0 4 — 127p 0 8 +72po 8 — 16p 0 10 )-E f 

n ** 2(1 -pa 1 ) 4 

_(3p 0 4 — 2po*+131p a 3 — 4p a 1<r )F' — (6po 2 — 31po 4 +108po , +53po 8 — 8po 10 )-E ,/ 

2/3X1 -Po 1 ) 4 


WU 


120/3 4 

C* 


WU 


TTl* 


(24po I0 4-609po , -462po e +273pa 4 -60p 0 !1 )i ^v 4-(120-561pa , + 1026po 4 -969p a # +48p a 8 -48p 0 l^l )£' ,, 

(1 — Po*) ! 

■jTTo _ P0 Z F' — PaE' 

" 0 - 0_ 2(1 — Pa 1 ) 

TO 

”°- 0_ (1-po 2 ) 

Tpu (3 po 8 — PtFjF' — (4 Po 8 — 2p 0 i )E' 

" 6(1 -po*)* 

Tpa 20* uro 

» i.o £5T H 2.J 


— „ (12p a *-7po 4 + 3 Po a )F'-(l 9p 0 *- 17 Po 4 + 6 P <?)E' 

24(1 — pa*) 3 

Tpa ^po^Spc^'— ^po^Po 4 )-^ 

24/3X1-P0 8 ) 4 

— a /3 2 [(7p 0 i +Po 4 )F , -(2+7p a *-po 4 )^'] 

U3 -°" 2(1 -pa 2 ) 3 

— , (p 9 4 +7p a 8 )F'+(po*-7po 4 -2po s )E' 

* 1_ 2(1 -pa 2 ) 8 

ii?. /3 4 [(5po , +3po 4 )F'-(l + 7p 0 2 )^'I 

Wx0 ~ (1-Po 4 ) 8 
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— « _ ^[(Spo’-Tpo^+^po^F'-ce-npoHiQpo 1 )^'] 

(1-P0 2 ) S " 

rn _ (60p 0 1 -44po 4 +44po e -12p 0 s )F'-(107p 0 ! -126po 4 +91po s -24p 0 s )£' 

120(1 -p,*)* 

Tfo _ (15p 0 4 +34p 0 8 -p 0 *)F'-(38p, 4 +12p 0 8 -2p 0 g )ff / 

81 1 20^(1 - Po 2 ) 4 

■ffa _ ^[ (33p fl t +14po 4 +p 0 c )F / -(6-h46po s -6po 4 +2p 0 8 )E1 
U a '° * *“ 


6(1 -Po 2 ) 4 

W* (6po 4 +44p 0 8 -2p 0 8 )F / +(3p 0 8 -34p 0 < -21p n 6 +4p o 8 )g / 

8,5 6(1 -po 1 ) 4 

WU=-WU 

ITi.,— in., 

W °*’* 

m _ 3^[( Po 3 + po 4 )F'-2 Po s ^'I 

M °- 0 2(1 — pn 4 )* 


IM 6ff 4 rrr 0 

1 0, 0 ✓TS ” 2,2 


jyo TITO 

} } 1,0 — 20~ n fl,a 

Tf?.o=-^plrs. 


in. 0 =-3Tn, 0 

Combinations of these functions useful in writing the potential in concise form are as follows: 

p I Ql t Rj_ 

1 P}+ M 2+ M* 

where 


P’1 — Aq-~A>i 

Qi=o 

P, = 0 

P 2 == i — Ai2 

Qt=Q 

#2=0 

A=2 (-rio - 2^4i-|-^lj) 

q*=-4? 

2?,= 0 

1! 

~' r 

Q Ax 

Qt 2 

/?4=0 

Pi— g (-rio — 3 -4.1 + 3 -4.2 — A*) 

Cs = — 2 

R s =0 

P # =i(3.4,-.4 6 ) 

Qt — — q (ri.i /1 3 ) 

R t = 0 

P 7 =4 (A l -2A,+A i ) 
£ 

Q,—f 

i? T =0 

r s =^(2A,-A t ) 


P,=0 
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P 9 — (-4 0 4 A 1 + 6 A 2 — 4Ai+4g) 

Q» 4 (-4-0 — 2A 2 -r-4.4) 

7? 

P l0= 24 (6-4. a — 4-4 j+.4. 7 ) 

$ 10= — 4 (-4-1 — 2A 3 -{-.4. e ) 

n 4_i 

itio==Yf 


O — " "^' s 

7 ? 4. 2 

11 24 

vn— — 4 - 


Pl 2 = B0 (Aa “ 5 Al+ 10 - 4l ~ 10 ^+5-4e- A s ) 

Qm ^9 (-4.0 3z4 2 -}-3j4i 4. 7 ) 

i2a=4 (-4. 0 — .4. s ) 

Pl 3 == r >0 ^0A[+oA 7 — 4. t0 ) 

n — l~(A~ qI 4.9 I X) 

jBi3=Y 

Vis Yo [ 0 VI 5 — -ris) 

Pu= v>0 (^^ s_ -^-u) 

Qu = — YW (3-4-8 — -4g) 

^u =4 (-4i — 4. s ) 

■f > i 6 = 0 ' (-4i — 3 Az + 3 Ai — rils) 

<2l5 = — 9 (^-2 — ^- 4 ) 

5]j=0 

Pts = ~Q (3.d. 5 — 34.* — A.?) 

Qlt = — 9 ^- 4.3 — 4. 5 ) 

b h =o 

-i 

11 

1 

®l k- 

Q 4-6 

Vit 2 

Pi-=o 

Pii—Tf^ (-4] — — 4.4b+.4 8 ) 

Qi»= — 4 (-4.S — 2-4.1+.4.;-) 

T? 4-3 

^*=- 4 - 

P w == ' 9 ^‘ C-^io — 44. 7 -{-64^ — 4.4.5) 

Qia = — y (■^•* — 24. s -t-4L s ) 

E> -4_1 

K ™=~± 

^ >20 = O 4 (^ u — 4,4*) 

Qm=l {2A t -A t ) 

n A.Q 
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T 


p * 

A % 

A , 

At 

0 

1.0000 


0.1667 

.03403 

.9074 

.0020 

.1030 

.08710 

.0876 

.0000 

.1682 

.17366 

.9616 

.0048 

.1811 

.25882 

.0200 

.8348 

.1100 

.84202 

.8041 

.7080 

.0023 

.42202 

.8603 

.7078 

. .0781 

.60000 

.8287 

.0662 


.67368 

.7046 

.0101 

.0680 

.04279 

.7000 

.6717 

.0610 

.70711 

.7404 

.6304 

.0466 

.76604 

.7178 

.6192 

.0411 

.81016 

.6082 

.4807 

.0370 


.6814 

.4712 

.0348 

.00631 

.6076 

.4602 

.0327 

.06603 

.0478 

.4366 

.0208 


.0360 

.4244 

.0283 


Pt p» 


Pt -Pi Pa Pit 


0 

-.00008 

-.00108 
.00016 
.00481 
.00981 
.01309 
. 01 Y 07 
. 0101 ? 
.09080 
.02149 
.02103 
.09220 
.09232 
.02238 
.02230 
.02223 


0 

-.000001 

-.000096 

.000227 

.000008 

.001076 

.001310 

.001400 

.001632 

.001664 

.001040 

.001640 

.001684 

.001601 

.001477 

.001430 

.001416 


0.003126 

.002088 

.002670 

.001878 

.001361 

.000087 

.000730 

.000668 

.000447 

.000368 

.000293 

.000261 

.000210 

.000101 

.000171 

.000147 

.000136 


.0000011 

.0000018 

-.0000084 

-.0000280 

-.0000410 

-.0000420 

-.0000207 

.0000007 

,0000270 

.0000678 


0 

-. 0000020 
-.0000178 
-.0000230 
.0000177 
.0000723 
.0001206 
.0001608 
.0001833 
.0002031 
.0002160 
.0002224 
.0002276 
.0002811 
.0002831 
.0002346 
.0002348 


.0000442 

.0001443 

.0002302 

.0002472 

.0002276 

.0001094 

.0001716 

.0001470 

.0001804 

.0001186 

.0000089 

.0000880 

.0000804 

.0000740 

.0000000 

.0000012 


-0.012600 

-.011061 


-.006408 

-.008040 

-.002066 

-.002278 


-.000702 

-.000086 

-.000680 

-.000630 


0 

-.0000064 

-.0000242 

0000738 

-.0000064 

-.0000600 

.0000376 

.0001283 

.0002002 

.0002770 

.0008288 

.0003716 

.0004000 

.0004231 

.0004388 

.0004676 

.0004064 


0 

-.000047 

-.000006 

.000342 

.000000 

.000710 

.000717 

.000070 

.000616 
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T l " 


Co) Q/ 


p * 

Qt 

Qi 

Qt 

Qt 

Or 

Ui 

Qt 

Qu 

Qn 

Qu 

Qu 

Qu 

Qu 

Qu 

Qll 

Qu 

Ql. 

Q» 

0 

0 

-0.08334 

0 

0 

0 

0.08334 

0 

0 

-0.00626 

0 

0 

0 

0 

0 

0.01260 

0 

0 

-0.006260 

.03403 

.00107 

-.08277 

,00011 

-.00007 

.00006 

,.08180 

.00001 

.00016 

-.00007 

.000002 

.000014 

-.000007 

.00009 

-.00120 

.01216 

.000007 

.000118 

-.006782 

.08710 

.00473 

-.08073 

.00078 

-.00408 

.00396 

.07660 

.00000 

.00040 

-.00661 

.000010 

.000060 

-.000348 

1 HUJ.ilil 

-.00488 

.01102 

.000068 

.000196 

-.004474 

.17300 

.01237 

-.07000 

.00320 

-.01001 

.00917 

.00694 

.00047 

.00011 

-.00444 

.000088 


-.000064 

.00220 

-.01023 

.00887 

.000316 

-. 000308 

-.002475 

.20882 

.01033 

-.07068 

.00647 

-.01492 

.01286 

.00674 

.00121 

-.00087 

-.00361 

.000161 

.000186 

-.000773 

.00406 

-.01018 

.00702 

.000724 

-.001272 

-.001189 

.34202 

.02474 

-.06840 

.00975 

—.01880 

.01400 

.04700 



—.00270 

.000333 

-.000016 

-.000782 

.00649 

-.01433 

.00568 

.001136 

-.001088 

-.000443 

.42202 

.02804 

-.06072 

.01280 

-.02007 

.01004 

.04000 

.00307 

-.00312 

-.00225 

.000630 


-.000742 

.00646 

-.01442 

.00440 

.001487 

-.002464 

-.000021 

.00000 

.03133 

-.00641 

.01490 

-.02203 

,01043 

.03438 

,00393 

-.00403 

-.00184 

.000740 

-.000431 


.00706 

-.01390 

.00368 


-.003740 


.87308 

.03312 

-.00209 

.01068 

-.02274 

.01644 

.02084 

.00466 

-.00474 

-.00168 

.000920 

-.000624 

-.000623 

.00739 

-.01328 

.00307 

.001880 

-.002869 

.000880 

.64270 

.03427 

-.04024 

.01802 

-.02304 

.01020 

,02820 

.00623 

—.00620 

-.00130 

.001086 

-.000788 

-.000660 

.00760 

-.01261 

.00260 

.001076 

-.002001 


.70711 

.03407 

-.04034 

.01000 

—.02304 

.01007 

.02320 

.00500 

-.00661 

-.00112 

.001217 

-.000923 

-.000622 

.00701 

-.01170 

.00224 


-.002873 

■ 


.03036 

-.04384 

.01071 

-.02288 

.01560 

.02007 

.00606 

-.00600 

-.00098 

.001324 

-.001011 

-.000486 

.00761 

-.01107 

.00100 


•Blu/ .1 


.81010 

.03600 

-.04173 

.02022 

-.02261 

.01633 

.01911 

.00033 

-.00500 

-.00087 

.001411 

-.001112 

-.000469 

.00767 

-.01047 

.00174 


-.002778 


.80003 

.03000 

-.08000 

.02068 

-.02231 

.01603 

.01766 

.00063 

-.00821 

-.00079 

.001477 

-.001166 

-.000410 


-.00080 

.00169 


-.002714 

.000446 

.00031 

•0300 D 

-.03801 

.02083 

-.02201 

.01476 

.01649 

.00069 

-.00028 

-.00078 

.001628 




-.00046 

.00146 


-.002061 

.000436 

.00090 

.08047 

-.03647 

.02110 

-.02102 

.01426 

.01610 


-.00636 

-.00006 

.001600 


-.000354 



.00130 

.002002 

-.002649 

.000416 


.03037 

-.03037 

.02122 

-.02122 

.01416 

.01406 

,00097 

-.00637 

-.00001 

.001628 

-.001201 

-.000337 

.00728 

-.00849 

.00121 


-.002492 

.000404 


< d ) Hi 

o . 

tit 


fin 

i?il 


fin 

fill 

iZja 

fi « 

0 

0 

0 

0.0081260 

0 

0 

0 

0 

0 

0.003126 

.08403 

.0000070 

-.0001188 

.0030020 

.0000027 

-.000014 

.0000663 

.0000044 

-.000106 

.003037 

.08716 

.0000400 

-.0004802 

.0029706 

.0060101 

-.000007 

.0002206 

.0000306 

*—.000802 

.002769 

.17366 

.0001060 

-.0011024 

.0027328 

.0000870 

-.000368 

.0004003 

.0001000 

-.000708 

.002233 

.26882 

.0003810 

-.0016841 

.0024780 

.0001031 

-.000070 

.0000880 

.0001870 

-.001008 

.001780 

.34202 

.0006620 

-.0020201 

.0022460 

.0003060 

-.000060 

.0008200 

.0002400 

-.001080 

.001426 

.42202 

.0000037 

-.0022367 

.0020406 

.0004100 

-.001161 

.0008864 

.0002831 

-.001074 

.001166 

.60000 

.0008010 

~. 0023441 

.0018022 

.0004060 

-.001813 

.0000128 

.0003047 

-.001031 

.000040 

.67368 

.0008814 

-.0023884 

.0017082 

.0006086 

-.001416 

.0009100 

.0003158 

-.000072 

.000791 

.64270 

.0000891 

-.0023026 

.0016768 

.0006188 

-.001481 

.0000075 

.0003206 

-.000012 

.000600 

.70711 

.0000799 

-.0023717 

.0014062 

.0006680 

-.001810 

.0008002 

.0003213 

-.000863 

.000876 

.76004 

.0010087 

-.0023897 

.0013724 

.0006884 

-.001638 

.0008002 

.0003203 

-.000802 

.000603 

.81015 

.0010278 

-.0022009 

.0012018 

.0007100 

-.001648 

.0008606 

.0003180 

-.000756 

.000441 

.86603 

.0010400 

-.0022660 

.0012272 

.0007203 

-.001848 

.0008200 

.000314$ 

-.000711 

.000408 

.00031 

.0010401 

-. 0022170 

. 0011742 

.0007370 

-.001640 

.0008046 

,0003113 

-.000678 

.000370 

.06603 

.0010670 

-.0021678 

,.0011006 

.0007633 

-.001630 

.0007747 

.0003041 

-.000030 

.000320 

1.00000 

.0010610 

-.0021221 

.0010610 

.0007579 

-.001610 

.0007670 

.0008032 

-.000006 

.000303 
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